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Radical-ion-pair reactions were recently shown to manifest a host of non-trivial quantum effects 
accounted for by quantum measurement theory. An alternative approach purporting to describe 
the fundamental quantum dynamics of spin-selective radical-ion pair reactions was introduced most 
recently, bringing to three the competing theories, including the one traditionally used in spin chem- 
istry. We here consider entropy as a fundamental concept enabling a comparison of the predictions 
of these theories against what is physically acceptable on quite general grounds. 
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Spin-selective radical-ion-pair reactions are at the core 
of spin chemistry 1]. It was recently shown Q that 
radical-ion-pair reactions manifest non-trivial quantum 
effects that were masked by the traditional treatment of 
these reactions used since the advent of spin chemistry. 
Using the specific Hamiltonian interactions pertaining in 
radical-ion pairs, specifically taking into account elec- 
tron tunnelling during charge recombination, the den- 
sity matrix equation describing non-reacting radical-ion 
pairs was derived |2j based on basic quantum measure- 
ment theory considerations. Most recently, yet another 
density matrix equation supposed to describe radical-ion- 
pair reactions was suggested by Jones & Hore |3|, advo- 
cated to also follow from quantum measurement theory. 
Moreover, two papers appeared [H, H| purporting to have 
derived from first principles the traditional master equa- 
tion of spin chemistry. We will here consider the purity 
of the spin state of radical-ion pairs as a guiding con- 
cept enabling a comparison of all three theories against 
what is expected on general grounds, like the second law 
of thermodynamics. We will not address the traditional 
theory individually, since it suffers from the exact same 
problems as the Jones-Hore theory, problems to be out- 
lined in the following. 

The fate of radical-ion-pair reactions is determined by 
the time evolution of /?, the density matrix describing the 
spin state of the molecule's two electrons and any number 
of existing magnetic nuclei. In 2009 Kominis derived the 
Lindblad-type master equation 
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where % is the magnetic interaction Hamiltonian induc- 
ing singlet-triplet mixing, ks and fcr are the singlet and 
triplet recombination rates and Qs and Qt are the sin- 
glet and triplet projection operators (only Qs appears 
in (p} because both the singlet and the triplet reservoir 
perform a continuous measurement of Qs)- The above 
trace-preserving master equation describes the density 
matrix evolution of non-reacted (i.e. not yet recom- 
bined) radical-ion pairs. It has been augmented by 
trace-decaying terms describing the recombination pro- 
cess @, but they are not needed for the present discus- 
sion. The master equation ([1} describes singlet-triplet 
decoherence induced by the continuous measurement of 
the singlet and triplet reservoirs. These reservoirs are the 



excited vibrational levels of the singlet and triplet neu- 
tral molecules resulting from the charge recombination, 
as will be illustrated in the following. 

A different approach was introduced by Jones & Hore 
[3j, based on the master equation 



(ks + kr)p + ksQrpQr + krQspQs (2) 



We will here examine the repercussions of this theory 
regarding the purity of the spin state of radical-ion pairs. 
For simplicity, we will here consider an idealized radical- 
ion pair with no nuclear spins and at zero magnetic field. 
That is, we completely neglect singlet-triplet mixing in- 
duced by magnetic interactions. Furthermore, we con- 
sider just one recombination channel, e.g. the singlet, 
and thus set kr = 0. Finally, we take as initial state the 
coherent superposition \ip) = (\S) + |T))/V2- So in this 
example the Jones-Hore master equation reads 



dp 

— = -k s p + ksQrpQr 
dt 



(3) 



The philosophy of postulating such a master equation 
is the following if the density matrix of a single 
radical- ion pair at time t is pt, then there are three 
different scenarios for what can happen in the follow- 
ing time interval dt: (i) with probability 1 — ksdt noth- 
ing happens, keeping pt+dt = Pt, (ii) with probability 
/c5<i<Tr{Qspt} a singlet recombination takes place leading 
to pt+dt — 0, and (iii) with probability /csG?iTr{QTPt} the 
molecule's state is projected to the triplet state, making 
Pt+dt = QTPtQT/Tr{Q T pt}- Averaging all three pos- 
sibilities leads to However, equation ([3]) was not 
derived from the usual procedure of considering (i) the 
quantum system under study and (ii) the measurement- 
reservoir degrees of freedom, evolving the combined state 
under the specific system-reservoir coupling and then 
tracing out the reservoir degrees of freedom. It was pos- 
tulated to follow from general considerations in analogy 
with the well known textbook example of a single two- 
level atom observed by a photodetector. If the atom is 
initially in a coherent state (|e) + \g))/y/2, and if the pho- 
todetector does not register a photon until time t, the 
state of the atom becomes (e~ r */ 2 |e) + \g))/Af, where 
M is the normalization constant and T the spontaneous 
decay rate. The interpretation is that not observing a 
photon, our knowledge of the atom's state acquires more 
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FIG. 1: Time evolution according to the Jones-Hore theory of 
(a) the singlet expectation value, (b) the triplet expectation 
value, (c) the number of existing radical- ion pairs, and (d) the 
purity of their quantum state, for the simple example of no 
singlet-triplet mixing, just the singlet recombination channel 
and a coherent initial state. 



ground-state " character" . The Jones-Hore theory makes 
the seemingly innocent analogy that in the case of a 
radical-ion pair that reacts through the singlet channel, 
the state of a non-reacting radical-ion pair acquires more 
triplet-state "character". However, this analogy breaks 
down if examined more carefully. Before explaining why 
this is the case, we will present the solution of ^ for the 
specific example we are considering. 

In FigQ] we plot as a function of time (a) the expec- 
tation value (Qs) = Tr{Qsp}, (b) the expectation value 
(Qt) = Tr{Qx/o}, (c) the normalization of the density 
matrix Tr{p}, which represents the number of surviving 
radical-ion pairs, and (d) the purity of the state Tr{p 2 }. 
The interpretation of these results is that (a) half of the 
radical-ion pairs react through the singlet channel, (b) 
the other half remain locked in the non-reacting triplet 
state, (c) the number of radical ion pairs is, as expected, 
reduced by half as f -> oo, and (d) we start with a pure 
state and end up with a pure (triplet) state. The lat- 
ter conclusion can be better visualized by considering 
non-reacting radical-ion pairs. To that end, we normal- 
ize (Qs) and (Qt) with Tr{p} whereas we normalize the 
state purity Ti{p 2 } with Tr{p} 2 . The results are shown 
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FIG. 2: Time evolution according to the Jones-Hore theory 
of (a) the singlet expectation value, (b) the triplet expecta- 
tion value, and (c) the state purity of non-reacting radical-ion 
pairs. 



in Fig|2] The interpretation is that as long as a radical- 
ion pair does not react, its state approaches the pure 
triplet. From the construction of the Jones-Hore master 
equation ([3]) , in the example at hand there are the follow- 
ing two possibilities for a non-reacting evolution between 
t and t + dt: (i) p t+dt = p t , or (ii) p t+dt = \T)(T\. The 
average of all possible trajectories embodying these two 
possibilities is what is plotted in FigJ5] So we start from 
a pure state, this becomes mixed, and then it returns to 
be pure again. At first glance this result is intuitive and 
acceptable. It seems to be analogous to a situation where 
a decaying atom is initially in the pure excited state, 
it becomes entangled with the radiation field, hence the 
atomic state is mixed, and then it decays to the ground 
state, becoming pure again. Or, if we have an ensemble 
of spins initially in the coherent state (| t) + I 4))/v2, 
and they undergo a Ti process, the state becomes mixed, 
| t>(t 1/2 + U)U I/ 2 - If a T i process follows, we would 
end up with a pure spin state again. However, here there 
is a problem. Among the N molecules, N/2 non-reacting 
ones will evolve in trajectories that average to FigJ2J So 
for the ensemble of non-reacting radical-ion pairs we have 
the entropy initially increasing (Tr{p } falls), and then 
decreasing (Tr{p 2 } rises back to 1). What makes up for 
this entropy decrease? The molecules are freely evolv- 
ing, i.e. there is no energy exchange nor any entropy 
exchange with the environment. These molecules do not 
react, so there is nothing to take up the entropy loss, 
as in the case of decaying atoms or spins. For example, 
when atoms spontaneously decay to the ground state, 
the entropy decrease mentioned above is counterbalanced 
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FIG. 4: Von Neumann entropy of non-reacting radical-ion 
pairs according to the Jones-Hore theory (blue solid line) 
and information gained by the reacting radical-ion pairs (red 
dashed line). The latter cannot accommodate the entropy 
loss taking place after t « 1/ks- 
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FIG. 3: Example of a quantum trajectories according to the 
Jones-Hore theory, (a) A reacting radical-ion pair according 
to the Jones-Hore theory. The expectation value of Qs jumps 
to one and at the same time ii the radical-ion pair singlet- 
recombines. (b) A non-reacting trajectory. At some random 
time ti the radical pair's state jumps to the triplet state and 
stays there forever. Blue solid line is the expectation value 
(Qs), red dashed line is Tr{p}. 



by the entropy increase of the radiation field. Similarly, 
when the spins undergo a T\ process and end up in a 
pure state after having been mixed by a Ti process, the 
entropy is taken up by the reservoir responsible for the 
T\ process. How is the loss of entropy of non-reacting 
radical-ion pairs balanced? 

To illustrate these considerations in more detail, we 
consider single quantum trajectories of reacting and non- 
reacting radical-ion pairs according to the Jones-Hore 
theory. As mentioned previously, the density matrix p t 
of a radical-ion pair at time t can evolve in either of the 
three ways (i) p t+dt = p t , (ii) Pt+dt = \T)(T\ and (hi) 
Pt+dt = 0. Let us first look at the evolution of a radical- 
ion pair that does react at some point, i.e. at some ran- 
dom time possibility (hi) is realized. Of course, until 
then p t = po, i-e. possibility (i) is realized. An example 
of such a trajectory is shown in FigJ3^. At time t\ the 
radical pair's state jumps to the singlet state and recom- 
bines at the same instant. An example of a non-reacting 
trajectory is shown in Fig|3L. At some random time 
*2, (Qs) jumps to (Qs) — ((Qt) — 1) and the state 
is locked there forever. So the average of non-reacting 
trajectories up to time t includes trajectories that have 
remained at the initial state po up to time t and tra- 
jectories that have been projected to the triplet state at 
some earlier time. That's where the decreasing purity of 
the non-reacting trajectories stems from. According to 
the Jones-Hore theory, the weight of the former decays 
as e~ kst , whereas the weight of the latter approaches 1/2 
as (1 — e~ kst )/2. That's why ast->oo the purity of the 
state of non-reacting trajectories approaches 1. But what 



after all compensates for the unphysical spontaneous en- 
tropy loss of the non-reacting radical-ion pairs? Is it the 
information gained from the reacting molecules? Not re- 
ally. In Fig|3] we plot (i) the von- Neumann entropy of 
p nr , the density matrix of the non-reacting-radical ion 
pairs, which is St = — Tr{p nr lnp nr }, and (ii) the infor- 
mation gained up to time t from the reacting molecules. 
This is easily calculated as follows: Let ps = Tr{Qsp n r} 
be the singlet probability at time t of a radical-ion pair 
just before it reacts. The information content of this 
molecule's state is sj = — psinps — (1 ~Ps) m (l — Ps)- 
When the radical-ion pair reacts we have ps = 1 , hence 
the information gain is si itself. There are ksdtTr{Qgp} 
molecules reacting in the time interval between t and 
t + dt, so the total information gain up to time t is 
It = J (iifcsTr{Qsp}si. From Fig|4]it is seen that from 
1 rs 1/ks an d on, the von-Neumann entropy of p nr drops 
from 0.4 to zero, but the information gain is just about 
0.1, i.e. it cannot accommodate for this entropy loss. 

We will now analyze the root of this shortcoming of 
the Jones-Hore theory. We will show where the analogy 
with the single atom mentioned above fails. In Fig|5]we 
pictorially depict the atom's two energy levels and the ra- 
diation reservoir coupled to the excited state. An atom 
initially in a coherent superposition of the ground and 
excited states has the possibility of decaying. A photon 
will thus occupy one of the states of the radiation reser- 
voir coupled to |e). If no photon is observed up to time 
t, one can claim that the atom's state has more \g) char- 
acter. There is no other possibility. In other words, if 
a single occupation leaves the atomic system (spanned 
by \g) and |e)) and goes to the radiation reservoir, then 
this event is definitely observable in the outside world as 
a photon. From the non-observation of a photon an un- 
ambiguous inference about the atom's state can be made. 
However, the situation is not the same with radical-ion 
pairs. In FigJS^, we show the energy levels of a radical- 
ion pair and in Figjfjb we simplify this picture for use 
in the following. Evidently, we here have two reservoirs. 
One is the reservoir of the excited vibrational states of 
the neutral molecule. There is a second one. Each of 
those excited states is coupled to a radiation reservoir 
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FIG. 5: (a) Two-level atom initially in a coherent state (|e) + 
\g))/y2. The radiation field reservoir is coupled to the excited 
state \e). (b) If the atom decays, a photon occupies one of 
the radiation reservoir states. This is definitely observable. If 
the atom does not decay, then the probability that it is in the 
ground state increases. 

(or another phonon reservoir, it does not matter). 

How is the event of the charge recombination transmit- 
ted to the outside world? It is by the emission of the pho- 
ton by the excited vibrational level and the creation of the 
singlet neutral product DA. If the radical-ion pair is at 
time t in the coherent superposition (IS") + \T))/^/2 or in 
any state for that matter, and no photon is detected, i.e. 
no reaction takes place within the following time interval 
dt, the only inference to be made is that the singlet reser- 
voir state was not occupied in this time interval. No in- 
ference can be made about the radical-ion pair spin state. 
It is exactly this point where the analogy with the atom's 
case fails. In other words, the basis of the Jones-Hore 
theory is the erroneous association of the quantum mea- 
surement going on internally in radical-ion pairs due to 
the internal spin-state-measurement dynamics with the 
recombination event itself. To clarify this statement, in 
FigEK we depict a radical-ion pair initially in the coher- 
ent state (IS") + \T))/^/2, i.e. the occupation resides in 
the quantum system spanned by IS") and \T). If the occu- 
pation leaves the system and goes to the singlet reservoir, 
this is not observable (FigJTp). What is observable is the 
next possible step (FigJTb) , i.e. the decay of the reser- 
voir state and the creation of a photon and the ground 
state DA. But this is just one possibility. There is a sec- 
ond possibility (Fig(7p). The occupation tunnels back to 
| S) and the spin state evolution continues on from there. 
This possibility is described by the second-order pertur- 
bation theory from which follows the master equation 
(P} derived in 0]. It is these virtual processes that lead 
to singlet-triplet decoherence introduced in Q. In other 
words, during a possible trajectory of a single radical-ion 
pair, one can have many instances where the continu- 
ous measurement induced by the singlet reservoir (and 
the triplet reservoir when it exists) results in Qs = 1 or 
Qs = 0. The Joncs-Hore theory wrongly associates the 
former possibility with a concomitant singlet recombina- 
tion. If the triplet recombination channel was also open, 
then the Jones-Hore theory would wrongly associate both 



D*A 



b 




DA 

FIG. 6: (a) Radical-ion pair quantum dynamics. A photon 
excites the singlet neutral precursor molecule DA into D*A, 
and a charge transfer creates the radical-ion pair. The excited 
vibrational levels (DA)* of the neutral DA molecule form the 
measurement reservoir, which (i) acts as a measurement de- 
vice for the radical pair's spin state, and (ii) acts as a sink 
of radical-ion pairs, i.e. in the event of recombination, the 
electron tunnels into a reservoir state and a fast spontaneous 
decay results in the ground state DA (which is the singlet 
product) and a photon emission. Similar for the triplet reser- 
voir, (b) A simplified version of (a) in the case of just a 
singlet recombination channel. The (orange) reservoir on the 
left of the singlet reservoir depicts the radiation reservoir (or 
some other phonon reservoir) to which the excited vibrational 
states of the singlet reservoir are coupled. 
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FIG. 7: Radical-ion pair (a) initially in the coherent singlet- 
triplet state, (b) The electron might tunnel to a singlet reser- 
voir state, from which will follow (c) either a photon emis- 
sion and the production of a singlet neutral product i.e. the 
molecule DA (possibility 1), or it might tunnel back (possi- 
bility 2) to the singlet radical-ion pair state (virtual process) 
commencing again the spin-state evolution. 
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FIG. 8: A possible quantum trajectory of a radical-ion pair 
starting out from the singlet state, including singlet-triplet 
mixing for illustrative purposes. The two red circles are the 
times where we have a measurement outcome of Qs — in 
the first and Qs = 1 in the second. Finally, at some random 
time (black circle) the radical-ion pair recombines, and from 
that point on Qs = Qt = 0. 



the singlet (Qs — 1) and the triplet (Qs = 0) measure- 
ment outcomes with the singlet and triplet recombination 
event. The actual charge recombination event is, how- 
ever, completely independent and can take place at any 
time. This is shown schematically for a more general ex- 
ample than the one considered here in Fig|8j The kinks 
produced in a smooth state evolution by the jumps to 
Qs = 1 or Qt = 1 are responsible for singlet-triplet de- 
coherence. That is, the average of all trajectories like the 
one in FigJS] without a recombination event reproduces 
the master equation ([TJ derived in Q. If we calculate 
the state purity according to this master equation, we 
find that it is a decreasing function of time. In the par- 
ticular example we are considering here, a possible quan- 
tum trajectory for a reacting radical-ion pair is shown in 
FigE^. The radical-pair's spin state jumps to the sin- 
glet at time t\ (measurement outcome Qs = 1), from 
which it recombines at time ti. There is nothing against 
t-2 > tx, and this fact is what the Jones- Hore theory fails 
to account for. An example of a non-reacting trajectory 
is shown in Figj9jD. At some random instant £3, the mea- 
surement outcome is Qs — 0, and the state jumps to the 
pure triplet, where it remains forever. We stress that up 
to time ti both trajectories are non-reacting trajectories. 
The average of all such pure non-reacting trajectories re- 
produces the singlet-triplet decoherence embodied in (flj , 
which results in a monotonically decreasing state purity, 
as shown in Figj9j:. In other words, in this particular ex- 
ample the infinite-time state of non-reacting trajectories 
is Poo = \\S)(S\ + \\T)(T\ as opposed to Poo = \T)(T\ of 
the Jones-Hore theory. In our approach, it is easy to see 
that the von-Neumann entropy of non-reacting trajecto- 
ries monotonically increases. Hence the problem facing 
the Jones-Hore theory does not exist. 

We can recapitulate the following consistent picture of 
what really happens. Up to some time t, non-reacting 
trajectories involve pure-state trajectories of molecules 
that have been projected to the triplet state, and pure- 
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FIG. 9: Quantum trajectory of (a) a reacting radical-ion in 
which a measurement outcome of Qs = 1 takes place at time 
ti and the recombination reaction at a later time t-z- This tra- 
jectory is non-reacting until time ti- There is nothing against 
t2 > t\, and this fact is what the Jones-Hore theory fails to 
account for. (b) a non-reacting in which a measurement out- 
come of Qs = takes place at time £3. Blue solid line is 
the expectation value (Qs), red dashed line is Tr{p}. (c) the 
average of all non-reacting trajectories produces the decaying 
purity resulting from (fjj) for the example under study. 



state trajectories of molecules that have been projected 
to the singlet state at some time earlier than t and have 
not yet recombined. Until the latter do recombine, we 
have a mixed state of increased entropy. After they do, 
we end up with a pure triplet state. This entropy loss is 
counter-balanced by the entropy increase in the final pho- 
ton/phonon reservoir. It is through the reacting molecules 
that entropy is transferred from the radical-ion-pair sys- 
tem to the environment. In contrast, the Jones-Hore 
theory ignores the presence of the singlet vibrational- 
reservoir states of the radical-ion pair. This reservoir 
is in between the quantum system (radical-ion pair) and 
the outside world (photon reservoir). What goes on in 
this singlet reservoir is unobservable. Hence any infor- 
mation or the lack thereof in the outside world cannot 
be associated with the quantum system's state. How- 
ever, the Jones-Hore theory does exactly that. Hence 
in this description, the non-reacting trajectories involve 
pure state-trajectories still in their initial state, and pure 
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state trajectories that have been projected to the triplet 
state based on the non-observation of a reaction. The 
latter is, as explained before, an ambiguous inference on 
the radical-pair's spins state, and that's why it results in 
an unphysical entropy loss of non-reacting trajectories. 
The difference between our consistent physical picture 
and the Jones- Hore theory is not just a philosophical cu- 
riosity. It has measurable consequences on the reaction 
time. If we define the latter to be the time at which 



the radical-ion pair population has decayed to a specified 
fraction of its original value, then the two approaches 
will predict different reaction times, exactly because the 
Jones-Hore theory makes the association "measurement 
outcome=recombination" . We will further analyze reac- 
tion times in another manuscript. 

I would like to acknowledge several fruitful discussions 
with Prof. Ulrich Steiner. 
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